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Abstract. Inhomogeneous cosmological models have had significant success in 
explaining cosmological observations without the need for dark energy. Generally, 
these models imply inhomogeneous matter distributions alter the observable 
relations that are taken for granted when assuming the Universe evolves according 
to the standard Friedmann equations. Moreover, it has recently been shown that 
both inhomogeneous matter and pressure distributions are required in both early 
and late stages of cosmological evolution. These associated pressure gradients are 
required in the early Universe to sufficiently describe void formation, whilst late- 
stage pressure gradients stop the appearance of anomalous singularities. In this 
paper we investigate the effect of pressure gradients on cosmological observations 
by deriving the luminosity distance — redshift relations in spherically symmetric, 
inhomogeneous spacetimes endowed with a perfect fluid. By applying this to a 
specific example for the energy density distribution and using various equations 
of state, we are able to explicitly show that pressure gradients may have a non- 
negligble effect on cosmological observations. In particular, we show that a non- 
zero pressure gradient can imply significantly different residual Hubble diagrams 
for z < 1 compared to when the pressure is ignored. This paper therefore 
highlights the need to properly consider pressure gradients when interpreting 
cosmological observations. 
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1. Introduction 



The idea that we live in a Universe dominated by dark energy is based on a 
model dependent interpretation of cosmological observations. A key assumption 
invoked in these models is that, throughout its evolution, the Universe is both 
homogenous and isotropic on all scales and is therefore adequately described 
by a Friedmann-Lemaitre-Robertson- Walker (FLRW) spacetime. However, there 
now exists a wealth of literature that relaxes the homogeneity assumption to 
explain key cosmological observations without the requirement of dark energy. In 
particular, numerous authors have shown that the dimming of supernovae type 



change in the predicted luminosity distance - redshift relation (IDabrowski & Hendrv . 


1998 


: IPascual-SanchezL 199J 


i Celerierl 2000t Tomita. 2000l 2001allbl: Ieuchietal.. 


2002 


: Godlowski et a .. 


2004 


: Alnes & Amarzeuiouil 20061 Chune & Romanol 


2006; 


Alnes & Amarzeuioui . 


2007 


: lEnavist & Mattssonl. 2007: 1 Alexander et all 


2007; 


Boleikol 120081 Garcia-Bellido & Haue:b0elleL 2008a! bl; Zibin et al. 


. 2008t Yoo et al.. 


20081: lEnavist. 20081 Boleiko & WvitheS. 2009k Celerier et al. 


. 20091). 


Such 



inhomogeneous cosmological models have since been used to describe much more 
than just the dimming of SN la data. Indeed the implementation of such models 
may suggest we live next to a Gpc scale void (for a short review see lEllisl . 
2008). Alternatively, it has recently been shown that inhomogeneous cosmologies 



utilizing the pressure-free Lemaitre-Tolman models can also reproduce the SN l a 



(Garcia-Bellido & Haugb0elle 



results by including a "giant hump" at the present time (jCelerier et all 120091) . 
Moreover, nu merous future pro b es ha v e been suggested based on Baryon Acoustic 
Oscillations <|Boleiko fe WvitheL 120091; iGarcia-Bellido fc Haugboelld . l2009h . SN la 
measure ments at higher re dshifts ( Clifton et all l2008h . time-drift of cosmological 
redshifts (lUzan et al. L 20081). spectral distort ions of the cosmic microwave background 
power spectrum ( Caldwell fc; Stebbins . 2008) the kinematic Sunyaev-Zel'dovich effect 



:ebbins 
2008bh 
r et all 



the shape of the distance - redshift relation 



Clarkson et al l. 120081 ; iFebruarv et alll2009h and even the cosmic neutrino background 
Jia fc Zhangl . 20081 ) which may definitively determine the correct model of the 



Universe. 

There are many different approaches to implement inhomogeneous cosmological 
models including spatial averaging, p ertubative meth ods and the use of exact solutions 
of the Einstein field equations (see ICelerierl . 120071 for a review of these methods). 
The exact solution method is considered to be the "...most straghtforward and 



devoid of theoretical pitfalls" ([Celerierl 120071 ) as they can be used to represent both 
strong and weak inhomogeneities. However, these models are difficult to implement 
due to the limited type of known solutions with either complicated symmetries 
or matter distributions. As such, a majority of work to date has utilized the 
Lemaitre-Tolman class of solutions which describe a s pherically symmetric spacetim e 
endowed with an inhomogeneous distribution of dust ( Lemaitre L 19331: Tolmanl 19341) . 
Generalizations i nclude the use of Stephani spacetimes ( Clarkson fc Barrett! Il999l 
iGodlowski et all [2004. IStelmach fc Jakackal l2006l IDabrowski et all |2007[ )7 Szafron 
solutions (Moffat, 2008) as well as variou s swiss - cheese model s using Lemaitre-Tolman 



solut i ons (for exampl e iBr ouzakis et all l2007t iMarra et all 120071: I Alexander et al 



20071 : IBrouzakis et all 120081 : iBiswas fc Notarl l2008t iKolb et all l2009h and Szekeres 
solutions (jBoleikol . I2009T ) embedded in FLRW spacetimes. Each of these models have 
been utilized because exact solutions describing the evolution of the spacetime is 
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known, however this comes at a cost as it implies some sense of physical realism 
is relaxed. For example, the Stephani class of solutions have inhomogeneous 
pressure distributions but the density is only a function of the temporal coordinate, 
whereas t he Lema i tre-To lman spacetimes have an inhomogeneous density but zero 
pressure. Boleikol ( 2006 ) however, has shown that the formation of observed voids 
( Hovle fc VogelevL 20041 ) mav be highly dependent on inhomogeneous radiation that 



is non-zero at the surface of last scattering and soon thereafter. This was deduced by 
studying the spherically symmetric Ei nstein field equations en dowed with a perfect 
fluid with a radiative equation of st ate. IBoleiko fc Laskvl l|2008 | ) further s t udied these 
models to allay previous concerns (jrlellabv fc Lakel . Il985t IBoleiko et al. . 120051 ) that 
anomalous singularities, known as shell crossing singularities, would occur throughout 
the evolution of these inhomogeneous spacetimes. In fact, whilst shell crossing 
singularities may still occur in Lemaitre-Tolman spacetimes, the incorporation of non- 
zero pressure gradients delays the onset of singularity formation until considerably 
after structure formation has taken place. This implies that press ure gradients also 
play an important role in the late term evolution of the Universe ( Boleiko fc Laskvl . 
20081). For revi e ws of inhomogeneous cosmolog ical models the reader is referred to 
Krasihskil (jl997l ): ICelerierl (120071 ): iHellabvj (|2009h . 

The aforermentioned cosmological models attempt to explain observations of the 
late term evolution of the Universe through inhomogeneities in the baryonic matter. 
However, there exists many alternative approaches which also offer promising insights. 
In particular, if dark energy exists and is not a cosmological constant but some kind 
of dynamical field or negative pressure fluid, then it must posses inhomogeneities 
due to gravitational interactions with both itself and with clumps of baryonic and 
dark matter. The effect of these inhomogeneities on structure f o rmation can be non- 



neglible fe.g. | Basilakosl 120031; iLinder fc Jenkins!. 20031; iLinderl. 20051; jKuhlen et al 
20051 lAbramo et all . l2007l 120081 ; iMota et all 12008 : lAvelino et all 120081) and p ossibl 
requi res full non- linear evolution for a complete description ( Abramo et al. . 1200. 
2009). Therefore, in these models the geometry of the spacetime is necessarily 



inhomogeneous, implying the non-linearities in the Einstein field equations will play 
a pivotal role, in part due to the non-zero pressure gradients of the inhomogeneous 
dark energy. 

Non-zero pressure gradients enter the Einstein field equations through the energy- 
momentum tensor, and hence effect the overall geometry of the spacetime. This 
therefore effects the way in which we interpret various cosmological observations by 
changing the overall luminosity distance - redshift relation. The purpose of this article 
is therefore to calculate the luminosity distance - redshift relation for spherically 
symmetric, perfect fluid spacetimes where the pressure and energy-density are allowed 
to vary both in the temporal and radial directiorJ§|. These models generalize those 
of Lemaitre-Tolman spacetimes, which can be recovered by setting the equation of 
state such that the pressure vanishes. Moreover, these spacetimes can be applied to 
dynamical dark energy models by choosing an equation of state where the pressure is 
negative. 

This paper is set out as follows: In section [2] we show the general form of 

§ Such spherically symmetric, perfect fluid spacetimes wer e first considered by iLemait rd JlQ3 3f) who 
even a llowed for anisotropic pressure distributions. Later, Podurets (1964a b) and Misner & Sharp 
( 1964) considered the isotropic pressure case, and the equations have since been accredited in the 
literature to Misner and Sharp. However, throughout the present article we refer to these spacetimes 
as "Lemaitre" models to give the appropriate credit. 
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the sphericaly symmetric, perfect fluid spacetimes, including relating the metric 
coefficients back to familiar cosmological entities such as the local Hubble rate and the 
critical densities. In section[3]we derive the overall redshift relation for our spherically 
symmetric perfect fluid with an observer situated at the coordinate origin (i.e. r = 0), 
also the luminosity distance in terms of the metric coefficients, hence arriving at 
the luminosity distance - redshift relation. In section |4] we show an example of 
the change in luminosity distance - redshift relation for a specific example of an 
inhomogeneous cosmological model with differing equations of state, and we provide 
concluding remarks in section [5] Throughout the paper we use coordinates in which 
the speed of light is unity unless otherwise explicitly stated. Moreover, greek indices 
range over ... 3. 



2. Spherically Symmetric Perfect Fluid Spacetimes 

The e volution of a spherically symmetric perfect fluid was first studied by iLemaitre) 
( 1933ft . The metric in comoving coordinates can be expressed as 

ds 2 = -a 2 dt 2 + dr 2 + R 2 dn 2 . (1) 

l + E 

Here, a(t, r) > is the lapse function, R(t, r) is the area-radius coordinate and 
E(t,r) > —1 is the curvature parameter. Moreover, the two-spheres are denoted 
by dfl 2 := dO 2 + sin 2 9d(j) 2 and a prime denotes differentiation with respect to the 
radial coordinate, r. 

When the pressure and density are functions of the temporal coordinate only, the 
lapse function can be set to unity without loss of generality, and 

R-ta(t)r and E^-kr 2 . (2) 

This is then the standard FLRW spacetime, where k = 0, ±1 is the curvature index 
and a(t) is the scale factor. 

A perfect fluid stress-energy tensor takes the form 

T^u = (p + P)u fl u fl + Pg^ v , (3) 
where u M is the four-velocity of the fluid, which is given in coordinate form by 

We express the Einstein field equations as — kT^ — Ag^ u , where A is 

the cosmological constant and k = 8ttG. The Gt component of the Einstein field 
equations is 

E _^A =n m 
2 {l + E) aR> ' ( ' 

where an overdot denotes partial derivative with respect to temporal coordinate, i.e. 
E = dE/dt, and a prime denotes partial derivative with respect to radial coordinate, 
i.e. a' = da/ dr. In the dust limit, the second term in this equation vanishes, implying 
E = E(r). The above equation is therefore an extra equation that must be satisfied 
from both the dust and the homogeneous cases. The remaining Einstein field equations 
can be expressed as 

R 2 — Ea 2 2RR'-E'a 2 RE 

up + A, (5) 



o 



2 R 2 a 2 R'R a 2 R{l + E) 



A2 „^ ^9 2RRu 2Raa'(l + E) 77r 1N , x 

R 2 + 2RR - Ea 2 A '- = -R 2 a 2 UP - A . 6 

a R 
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It is again noted that, in the dust case, the final term on the left hand side of ([S]), 
and the last two terms on the left h and side of (Et vanish, and the equations reduce 



to those of the dust version given bv lEnqvistl H008), and indeed without these terms 



the equations are remarkably similar to the set of equations for homogeneous FLRW 
spacetimcs. Both the final terms on the left hand sides of equations §5§ and © are 
associated with the pressure gradient through equations (| 10[) and ((4]) . 
The first integral of equation ([6]) can be found to be 

R 2 F A E 



a 2 R 2 R 3 3 R 2 ' 

where F(t, r) is twice the Misner-Sharp mass (M isner Sz Shard . IT964t ) . This equation 



is now expressed analo gously to that of the FLRW and pressure-free Lemaitre-Tolman 



cases (for example see lEnqvist fc Mattssonl . 120071) . Indeed, the only differences, aside 



from the additional functional dependencies, is the presence of the lapse function, a, 
which we will later show becomes incorporated into the definition of the local Hubble 
rate. 

Putting this back through equations ([6]) and ([5]) respectively gives 

F =-kPRR 2 , (8) 

F' = k P R'R 2 . (9) 

When F' remains non-zero and R' or R tend to zero, equation ^ implies the 
density diverges, which is interpreted as shell crossing and shell focussing singularities 
res pectively (for a discussi on of shell crossing singularities in a cosmological context 



see 



Boleiko fc Laskv . 20081 and references therein). 



The radial component of the conservation of energy momentum, T a Ka = 0, 
implies 

(p + P)- = -P'. (10) 

a 

This equation says that the pressure gradient being non-zero, i.e. P' ^ 0, implies 
that the observer is non-geodesic. That is, the presence of a pressure gradient pushes 
the observer off the freely falling geodesies. If there is zero pressure gradient, then 
the lapse function is also a function of only the temporal coordinate, which can then 
be rescaled such that a — 1. This implies that an observer in a Universe with zero 
pressure gradient is freely falling along with the fluid. Moreover, equations (Q| and 
(|10l) imply that the temporal evolution of the curvature parameter, E, is governed by 
the pressure gradient. That is, if the pressure distribution were homogeneous (even if 
it were non-zero), then the curvature parameter would not evolve in time. 

The temporal component of the conservation equations can be manipulated using 
equation (fTU|) and (0| to give 

This equation has been expressed such that the explicit dependence on the pressure 
gradient can again be seen. In the Lemaitre-Tolman case, i.e. where P = 0, and the 
FLRW case, i.e. where P' — 0, the second term on the left hand side of the above 
equation vanishes. Moreover, in the FLRW case, the right hand side can be shown to 
be -3 (p + P) R/R. 
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We now define the local Hubble rate to be 



H(t, 



R 

aR' 



(12) 



The additional dependence on the lapse function, a, in this definition comes from 
equation (J7J , and is attributed to the observer no longer traveling along geodesies due 
to the introduction of a pressure gradient. Therefore, different observer will measure 
a different local Hubble expansion rat^jj. 

We define the current local Hubble rate, Hq(t) := H(to,r), defined on some 
spacelike hypersurface, given by t = to, and also the current area-radius coordinate, 
Ro(r) := R(to,r). In FLRW cosmology (with non-zero pressure contribution), the 
Hubble rate is given in terms of the critical densities according to 

? X -3(1+W) , jy X -2 



H 2 



Hi 



(13) 



where f2 c := 1 — il\ — (Im and tildes are used to denote quantities in the FLRW model. 
Analogously, from equation (UJ we make the following definitions for the various critical 
densities in terms of the metric coefficients in the Lemaitre models: 

F 



H l R l 
A 



H Q R 



1 - A - ^ 



This implies equation ([7J is now expressed as 
H 2 (t,r)=H*(r) 



H M I ^ 



R 
Rq 



-3(1+/) 



Kq 



-2(l+e) 



(14) 
(15) 

(16) 

^aW(17) 



where 



v ito/ Fo \Ro J Eq 

Thus, in analogy to FLRW spacetimes we see three qualitative differences. Firstly, 
we see that the fi's are no longer constant parameters but they depend on the radial 
coordinate, r. Secondly, / in (fTTf is no longer the equation of state as Cj is in (fT3"|) . 
Lastly, and most importantly, the function E depends on time (in FLRW models 
E/r 2 is constant), it is no longer constant and its evolution affects the expansion 
rate. This is an important observation since, for example evolving homogeneous dark 
energy with a given equation of state, w(z), will generally behave differently than an 
inhomogeneous dark energy fluid. 

In the Lemaitre- Tolman dust case equation (|13[) is separabl e and o ne can therefore 



-2e 



(18) 



find an implicit solution for R (see equation 2.22 of lEnqvistl . 12008 *). However, this 
is not the case when there exists a pressure gradient, implying we require numerical 
rather than analytic solutions. 

|| This is also an underlying pr e mise in the Fractal B ubble Cosmo l ogical models of Wiltshire and 
collaborators iWiltsh ire. 2007b a, 2008), although see iKwan et al.l ||2009|) for a recent critique of 
these models. 
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3. Distance Redshift Relations 

Herein we consider an observer at the origin of the coordinate system. The spherically 
symmetric nature of the spacetime implies all observed photons coming from our past 
null cone have come along radial null geodesies. Therefore, setting dO = dtp = ds 2 = 
in the line element ([T]) implies 

dt dr R' dr , . 

= •= yU r ) (19) 

du du oVl + E ' du ' n v ' 

where u parametrizes the geodesies and x nas been defined for simplicity in the 

following expressions. Now we consider two light rays separated by a small distance 

A(u), such that the first and second light rays are expressed as a parametrization of u 

as t±(u) = t{u) and ^(w) = t(u) + X(u) respectively. Differentiating these with respect 

to the parametrization, and substituting equation (j!9[) gives 

dti dr , . . . 

-T = —rX{t,r), (20) 
du du 

dt2 dr , . dX . . 

-T = -^-x(t,r ) + —. (21) 
du du du 

The equation for ti can also be calculated using a Taylor series approximation such 
that 

dt 2 dr 

= -^r.x(t + x,r) 



( , s . dx (t, r) ' 



(22) 



du du' 
dr 
du 

where terms non-linear in A have been neglected. This is a valid assumption as the 
separation of the two light rays is assumed to be much less than the distance over 
which we are considering the propagation of the light rays. Evaluating the right hand 
side of equations (f2"TT) and (|2"2")l gives a differential equation for A(tt), 

du du dt 
Finally, the redshift is defined as 

Differentiating this, and using the above relations implies 
dz dr dx ^ ^ 
du du dt 

_ dr (1 + z)R' ( R' a E N 

~~ du ay/1 + E yR' a 2^/1 + E / 

Finally, utilizing equation (fl9|) and (|25|) , one can derive a pair of differential equations 
which give the coordinates in terms of the observable redshiftfl 

dz x[l + z ) 
dz x(l + z) 

Here we have implicitly assumed that the redshift is a single-valued function of the geodesic 
parametrization, u. In Lemaitre-Tolman dust spaceti mes this has been show n to not necessarily be 
the case depending on initial density configurations tMusta pha et al. 

L Il998f) . and we note that it 
requires verification for the Lemaitre spacetimes with non-zero pressure distributions. 



(25) 
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where 

R' fit' a E \ 

X ~ aVTT^ \B/~a~ 2VTTE J ' ( ' 

The E term is related through equation (fT6"]) to the temporal evolution of the radiation 
and matter densities. That is, the way one views the redshift of an object depends on 
the evolutionary history of terms such as £Ir and VLm- 

We now derive relations for the distance measurements in terms of the spacetime 
variables. The area distance, Da, is found as the root of the ratio of the angle 
subtended by null geodesies diverging from the observer to the cross-sectional area of 
this bundle. For a general metric in spherical coordinates, null geo descis propagating 



along the constant 9 and <f> direction, the angular distance is (jEllis et all 1198 



Stoeger et all 119921 ) 



V 2 _ V|detj ab | 



Da= V ' ■ =R'(fr), (29) 

sm u 

where j a {, is the 9 — <p part of the metric. The luminosity distan ce, Dl, and the 
area distance are related by Di 2 = Da 2 (1 + z) A ( Etheringtonl . 19331 ) implying in the 
spherically symmetric perfect fluid spacetime, 

D L = R (1 + zf . (30) 



Equations ([26]). ([27]) and ([30]) now provide the necessary equations to determine the 
luminosity distance - redshift relation for specific matter distributions in a Lemaitre 
spacetime with a central observer. 

Finally, the luminosity distance allows us to work out the apparent bolometric 
magnitude, to, of a standard candle from the absolute bolometric magnitude, M, at 
a given redshift. This is expressed in Mpc as 

m = M + 5 logAc (z) + 25 (31) 

where Dl depends on all f2 a terms and the local Hubble rate. This expression is 
required in section 0] below for calculating the residual Hubble diagrams. 



4. Inhomogeneous Example 

We now provide an example of a calculation for an inhomogeneous cosmological model. 
To determine the effect of the pressure gradient on the luminosity distance - redshift 
relation, we concentrate on two models: the pressure-free Lemaitre- Tolman model and 
the Lemaitre model that we have been dealing with hitherto. We specify both models 
with the same initial data at the current insant. Moreover, the background model is 
chosen to be an open Friedmann model with Q ma t — 0.3, fl\ = and Hq = 67.2 km 
s _1 Mpc" 1 . The radial coordinate is defined as the areal distance, R, at the current 
instant: f := R(to,r). However, for clarity in further use, the tilde sign is omitted 
and the new radial coordinate is referred to as r. The pressure-free Lemaitre- Tolman 
model is of a Gpc-scale void model, with Eq and Fo [note from equations (jTUJ) , (0]), 
and ([S]) that for pressureless matter these functions depend on r only] calculated in 
the following way. The Fq function is found from ([9]) with p given by: 

p(t ,r) = p b 



1 + S — S exp 



(32) 
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where pi, = 0.3 x (3Hq)/ (8irG), S = 1.5, and a = 1 Gpc. The function Eq is evaluated 
by integrating (J7|) and assuming the the cosmological and integration constants are 
zero: 

R 

I ; ^ = ct . (33) 

J ^E + F /r 

Initial data for the Lemaitre model is E(to,r) = Eg and F(to, r) = Fq. In addition we 
assumed the following equation of state: 

P = Kp n , (34) 

where K and n are constant. We then use eqations ([26 j) - (|27]) to find t and r along the 
past null cone and simultaneously solve (J4j) , ([Tj) - and ([8]) to find the evolution of the 
system. 

Our results are presented in figures HHH Figures [1] and [5] present the calculated 
angular diameter distance - redshift and residual Hubble diagram respectively for 
different values of the equation of state parameter K . Figures [3] and [4] are respectively 
the same figures where we now vary the parameter n. 

It is apparent that for a linear equation of state, an increase in the proportionality 
constant, K , results in increasingly larger differences in the distance measurements 
(figures Q] and [5]). For low redshifts, z < 1, the pressure- free Lemaitre- Tolman model 
has smaller values of the distance at equal redshifts. Moreover as K increases, so 
too does the value of the distance. This feature is better depicted in the residual 
Hubble diagram, figure [21 where the variation in magnitude as a function of redshift 
is large at small redshifts but becomes less distinguishable for z > 1. Deviations 
from a linear equation of state are depicted in figures [3] and SJ In these figures we 
have only varied n between 1.05 and 0.95, thereby indicating the sensitive nature of 
this parameter on the distance measurements. This is due to small variations in n 
implying large variations in the magnitude of the pressure gradient with respect to 
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Figure 2. Residual Hubble diagram for models with equation of state P = Kp. 
The dashed curve has K = 0, and the solid curves from the top to the bottom (at 
the maximum) have K = 0.005, 0.002, 0.001 and 0.0005. 



o 

Q. 

< 

Q 




Figure 3. Angular diameter distance as a function of redshift for models with 
equation of state P = 10 
have n = 1.05, 1.00, 0.97 and 0.95 



3 p n . The curves from top to bottom (at the maximum) 



the energy-density. Again, the residual Hubble diagram highlights the importance of 
this parameter at low redshifts, z < 1, where a variation of n by 0.1 can cause Am to 
vary by up to 0.6. 

Another interesting feature is that the position of the maximum of the angular 
diameter distance shifts to low er redshifts as the amplit ude of the pressure gradient 
increases. As recently shown bv lAlfedeel Sz Hellabv ( 2009h the position of the apparent 
horizon (d-D/i/dz = 0) for non-zero pressure models is given by the same relation as 
for the dust models, i.e. R — 2M. However, since for p' ^ the function M depends 
on time and thus this relation is satisfied by a different pair of t and r than in the 
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Figure 4. Residual Hubble diagram for models with equation of state P 



10 _J p™. The curves from top to bottom (at the maximum) have n = 0.95, 
1.00 and 1.05. 



0.97, 



dust models. 

As can be seen from figures HHH the incorporation of pressure gradients in 
cosmological models can significantly alter the distance - redshift relations. For models 
with a positive pressure gradient, the luminosity distance for a given redshift is larger 
compared with models for P = 0. It is therefore apparent that cosmological models 
with non-trivial pressure distributions will yield different results than cosmological 
models that do not consider pressure gradients. 

5. Conclusion 



Whilst a 
feet of i 



considerable amount of 



literature 
il models 



has 
to 



been dedicated to the ef- 
explain the observed dis- 



tance redshift relations from SN la data (IDabrowski & Hendrvl. 1998 


Pascual-Sanchezl 1999 


Celerier, 2000t Tomital. |2000L l2001allbl: Iguchi et al. 


20021 Godlowski et al.. 


2004 


: Alnes & Amarzguioui, 20061 Chung & Romanol 2006 


Aines & Amarzguiouil 


2007 


; lEnavist & Mattssonl. l2007t lAlexander et all 12007 


Boleikol |2008[. jGarcia-Bellido & HauEbtfelleL |2008allbl: IZibin et all l2008t lYoo et al. 


20081: EnavistJ. [20081: Boleiko & Wvithe. 2009: Celerier et al.. 20091). surprisingly lit- 



tie has been done on the effect of pressure gradients in such models. The most simple 
invocation of pressure gradients into inhomogeneous cosmologies is to utilise spheri- 
ca lly symmetric Lemaitre models with an observer placed at the center of symmetry. 
In Boleiko & Laskv (2008), we studied these models to show th at anomalous singular- 
ities that may occur in the evolution of such dust spacetimes ( Hellabv fc Lake . 1985t 



Boleiko et all 120051 ) are delayed by the inclusion of realistic pressure gradients until 
significantly after structure formation has taken place. In the present article we have 
continued this line of research to derive the distance - redshift relations in spherically 
symmetric spacetimes endowed with a perfect fluid. Moreover, we worked in comov- 
ing coordinates with the observer at the center of symmetry. This work will therefore 
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allow for an interpretation of cosmological data using such models. 

By exploiting specific examples of inhomogeneous matter distributions with 
various "simple" equations of state, we have shown that the introduction of pressure 
gradients can have a non-negligble effect on observations, particularly for z < 1 
(figures [TH3]). This implies that one should definitely consider the full effect of pressure 
gradients in interpreting cosmological data-sets. However, this is a non-trivial task 
due to the increased number of free-parameters in the system. Indeed, if the data can 
be realistically fit using a void model with P = 0, then it can equally well be fit by 
implementing slightly different matter distributions with P ^ 0. The problem then 
becomes one of reverse engineering the correct spacetime geome try from the numerou s 
observations, a task which has been discussed in great detail in iLu fc Hellabvl (l2007h : 
McClure fc Hellabvl l|2008h . 



We conclude by remarking that this work may also have a significant effect on 
interpretations of models with dark energy which is not a cosm ological co n stant , 



but instead a dynamical field or negative pressure fluid (e.g. Basilakosl | 200 



Linder & Jenkins, 2003: Linder, 
Mota et al.l . l2008HAvelino et all 



2005 



2008; 



Kuhlen et ail 120051: lAb ramo et all 12007. 200 



Abramo et al.L [2009h . These models necessarily 



have inhomogeneous spacetime geometries, which is equally attributed to "pressure" 
gradients of the dark energy field. In this way, the distance - redshift relations that we 
have derived herein will be equally valid, and can therefore also be used to interpret 
cosmological data within such models. 
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